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Abstract 

We numerically work out the impact of the general relativistic Lense-Thirring effect 
on the Earth- Mercury range \p\ caused by the gravitomagnetic field of the rotat- 
ing Sun. The peak-to peak nominal amplitude of the resulting time-varying signal 
amounts to 1.75 x 10 1 m over a temporal interval At = 2 yr. Future interplanetary 
laser ranging facilities should reach a cm-level in ranging to Mercury over compa- 
rable timescales; for example, the BepiColombo mission, to be launched in 2014, 
should reach a 4.5 — 10 cm level over 1 — 8 yr. We looked also at other Newtonian 
(solar quadrupole mass moment, ring of the minor asteroids, Ceres, Pallas, Vesta, 
Trans-Neptunian Objects) and post-Newtonian (gravitoelectric Schwarzschild solar 
field) dynamical effects on the Earth-Mercury range. They act as sources of system- 
atic errors for the Lense-Thirring signal which, in turn, if not properly modeled, may 
bias the recovery of some key parameters of such other dynamical features of motion. 
Their nominal peak-to-peak amplitudes are as large as 4 x 10 5 m (Schwarzschild), 
3 x 10 2 m (Sun's quadrupole), 8 x 10 1 m (Ceres, Pallas, Vesta), 4 m (ring of mi- 
nor asteroids), 8 x 10 _1 m (Trans-Neptunian Objects). Their temporal patterns are 
different with respect to that of the gravitomagnetic signal. 
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1 Introduction 



In its slow-motion and weak-field approximation, the Einsteinian General 
Theory of Relativity (GTR) predicts that a slowly rotating central body of 

— * 

mass M and proper angular momentum S induces two kinds of small pertur- 
bations on the otherwise Keplerian motion of a test particle orbiting it. The 
largest one is dubbed "gravitoelectric" (GE) (Mashhoon, 2007), and depends 
only on the mass M of the body which acts as source of the gravitational 
field. It is responsible of the well-known anomalous secular precession of the 
perihelion u of Mercury of 42.98 arcsec cty -1 in the field of the Sun. There is 
also a smaller perturbation, known as "gravitomagnetic" (GM) (Mashhoon, 
2007), which depends on S: it causes the Lense-Thirring (LT) precessions of 
the node Q and pericenter a; of a test particle (Lense & Thirring, 1918) . More 
specifically, the Parameterized Post-Newtonian (PPN) perturbing acceleration 
App N to be added to the Newtonian monopole v4 Newt0 n = —GMf/r 2 in the 
equations of motion 

—j^ = ^Newton + ^PPN (1) 



is (Soffel, 1989, p. 89, p.95), (McCarthy & Petit, 2004, p. 106) 



A PPN = -E g - 2 x B g . (2) 



In it (Soffel, 1989, p. 89, p.95), (McCarthy & Petit, 2004, p. 106), 
E g = -§f { [m&M. _ V] r + 2(1 + 7 ) (r • v) v] }, 

(3) 

B 9 = -(^)£[SS(S-r)r], 



where 7, (3 are the usual PPN parameters (Will, 1993) which are both 1 in 
GTR, c denotes the speed of light in vacuum, G is the Newtonian constant of 
gravitation, and v is the velocity of the test particle moving at distance r from 
M; the unit vector f is directed from the central body to the test particle. In 
eq. (2)-eq. (3) E g is the GE field, while B g is the GM one. The resulting GM 



2 



precessions are, in GTR, (Lense & Thirring, 1918) 
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(4) 



where a, e, J are the semi-major axis, the eccentricity and the inclination to 
the equatorial plane of the central body, respectively, of the orbit of the test 
particle. The PPN expressions of the GM node and pericenter precessions have 
a multiplicative (1 + 7)/2 factor in front of (Ilt and wlt (Peterson, 1997, p. 
40). 

In regard to the LT effect, at present, some attempts to measure it in 
the gravitational fields of the Earth and Mars have been performed with the 
LAGEOS (Ciufolini et al., 2009) and Mars Global Surveyor (Iorio, 2006) ar- 
tificial satellites in a series of tests exploiting non-dedicated spacecraft in an 
"opportunistic" approach. Their status is somewhat uncertain, and the real- 
istic evaluation of the accuracy reached in such tests is matter of controversy 
(Krogh, 2007; Iorio, 2009; Ciufolini et al., 2009; Iorio, 2010) The data analysis 
(Everitt et al, 2009) of the Gravity Probe B (GP-B) mission (Everitt et al., 
2001), aimed to directly measure another GM effect, i.e. the Schiff (1960) pre- 
cession of the spin of an orbiting gyroscope, in a dedicated spacecraft-based 
experiment orbiting the Earth, has recently released its final results (Everitt 
et al., 2011); it was proposed for the first time in 1961 (Fairbank & Schiff, 
1961). Anyway, it will likely not be possible to repeat such an experiment in 
any foreseeable future because of its extreme sophistication and cost. 

Concerning the Sun and its planets, it can be said that, contrary to the 
planet-spacecraft scenarios, the systematic errors caused by competing clas- 
sical forces are of relatively less concern with respect to the measurement 
errors which, at present, represent the major obstacle in measuring the solar 
LT effect. Indeed, the non-gravitational perturbations are absent, while the 
gravitational ones are relatively well known and easily to model. Suffices it 
to say that, in contrast with, e.g., the Earth, only one even zonal harmonic 
coefficient J2 of the multipolar expansion of the non-spherical part of its New- 
tonian gravitational potential has to be taken into account for the Sun. For 
a long time since its prediction in 1918 within GTR, the LT precessions have 
been retained too small to be detected with planetary motions. Indeed, by 
assuming a homogeneous and uniformly rotating Sun, de Sitter (1916) ob- 
tained a value of —0.01 arcsec cty -1 for Mercury in his pioneering work in 
which he preliminarily worked out the effects of the solar rotation on the 
planetary perihelia within GTR limiting to ecliptic orbits only; such a figure 
is also quoted by Soffel (1989, p. 111). With the same assumptions concerning 
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the rotation of the Sun, Cugusi and Proverbio (1978) obtained —0.02 arsec 
cty -1 ; as a consequence, all such authors concluded that, at their time, it was 
not possible to measure the solar LT effect; Soffel (1989, p. 23), e.g., reports 
failed attempts to detect the LT effect by fitting planetary data with modern 
numerically produced PPN ephemerides. Nowadays, the expected magnitude 
of the LT planetary precessions is even smaller than before. It is so because 
recent measurements of the Sun's proper angular momentum 

S Q = (190.0 ± 1.5) x 10 39 kg m 2 s- 1 (5) 



from helioseismology (Pijpers, 1998, 2003), accurate to 0.8%, yield a value 
about one order of magnitude smaller than that obtained by assuming a ho- 
mogeneous and uniformly rotating Sun. Despite such new findings, the per- 
spectives of measuring the GM field of the Sun from planetary motions are 
more favorable now than in the past. This can be easily recognized from the 
following simple arguments. The characteristic length with which the accuracy 
of the determination of the orbits of the particles should be compared is 

^1^319*; (6) 



in the case of the GE effects, l g is usually replaced by the Schwarzschild radius 
R g = 2GM/c 2 = 3 km for the Sun. The present-day accuracy in knowing, e.g., 
the inner planets' mean orbital radius 



e 2 * 



{r) = all + -\, (7) 



is shown in Table 1. Such values have been obtained by linearly propagat- 
ing the formal, statistical errors in a and e according to Table 3 of Pitjeva 
(2008); even by re-scaling them by a a factor of, say, 2 — 5, the GM effects 
due to the Sun's rotation fall, in principle, within the measurability domain. 
Another possible way to evaluate the present-day uncertainty in the plane- 
tary orbital motions consists of looking at different ephemerides of compa- 
rable accuracy. In Table 1 we do that for the EPM2006/EPM2008 (Pitjeva, 
2008, 2010), and the DE414/DE421 (Standish, 2006; Folkner et al., 2008) 
ephemerides; although, larger than 5 (r) , the maximum differences between 
such ephemerides are smaller than the solar GM length Z®. 

Here we will focus on the Earth-planet range |pj because it is a direct, 
unambiguous observable for which great improvements are expected in the 
near future. 

The paper is organized as follows. In Section 2 we deal with the numerical 
calculation of the impact of several dynamical perturbations on the Earth- 
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Mercury range after having reviewed the extent of the expected future im- 
provements in measuring it. In particular, Section 2.1 and Section 2.2 are 
devoted to the general relativistic gravitoelectromagnetic fields of the Sun. In 
Section 2.3-Section 2.5 the effects of the Newtonian perturbations due to the 
solar oblateness, the ring of the minor asteroids, Ceres, Pallas and Vesta, and 
the TNOs are investigated. The conclusions are in Section 3. 



2 The interplanetary ranging 

Recent years have seen increasing efforts towards the implementation of the 
Planetary Laser Ranging (PLR) technique accurate to cm-level (Chandler et 
al., 2005; Smith et al., 2006; Neumann et al., 2006; Degnan, 2006; Turyshev & 
Williams, 2007; Merkowitz et al, 2007; Degnan, 2008; Zuber & Smith, 2008). 
It would allow to reach major improvements in three related fields: Solar Sys- 
tem dynamics, tests of GTR and alternative theories of gravity, and physical 
properties of the target planet itself. In principle, any Solar System body en- 
dowed with a solid surface and a transparent atmosphere would be a suitable 
platform for a PLR system, but some targets are more accessible than others. 
Major efforts have been practically devoted so far to Mercury (Smith et al., 
2006) and Mars (Chandler et al, 2005; Turyshev & Williams, 2007), although 
simulations reaching 93 au or more have been undertaken as well (Degnan, 
2006, 2008). In 2005 two interplanetary laser transponder experiments were 
successfully demonstrated by the Goddard Geophysical Astronomical Obser- 
vatory (GGAO). The first utilized the non-optimized Mercury Laser Altimeter 
(MLA) on the Messenger spacecraft (Smith et al., 2006; Neumann et al., 2006), 
obtaining a formal error in the laser range solution of 0.2 m, or one part in 
10 11 . The second utilized the Mars Orbiting Laser Altimeter (MOLA) on the 
Mars Global Surveyor spacecraft (Abshire et al., 2006; Neumann et al., 2006). 
A precise measure of the Earth-Mars distance, measured between their cen- 
ters of mass and taken over an extended period (five years or more), would 
support, among other things, a better determination of several parameters of 
the Solar System. Sensitivity analyses point towards measurement uncertain- 
ties between 1 mm and 100 mm (Chandler et al., 2005). Concerning Mercury, 
a recent analysis on the future BepiColombo 2 mission, aimed to accurately 
determining, among other things, several key parameters of post-Newtonian 
gravity and the solar quadrupole moment from Earth-Mercury distance data 
collected with a multi-frequency radio link (Milani et al., 2002, 2010), points 
toward a maximum uncertainty of 4.5 — 10 cm in determining the Earth- 
Mercury range over a multi-year time span (1-8 yr) (Milani et al., 2002; Ashby 

2 It is an ESA mission, including two spacecraft, one of which provided by Japan, 
to be put into orbit around Mercury. The launch is scheduled for 2014. The con- 
struction of the instruments is currently ongoing. 



5 



et al., 2007; Milani et al., 2010). A proposed spacecraft-based mission aimed 
to accurately measure also the GM field of the Sun and its J2 along with 
other PPN parameters like 7 and f3 by means of interplanetary ranging is the 
Astrodynamical Space Test of Relativity using Optical Devices 3 (ASTROD) 
(Ni, 2008). Another space-based mission proposed to accurately test several 
aspects of the gravitational interaction via interplanetary laser ranging is the 
Laser Astrometric Test of Relativity (LATOR) (Turyshev et al, 2009); the 
GM field of the Sun is one of its goals. 

We will concentrate on Mercury both because the magnitude of the GM 
signal is the largest one and in view of the aforementioned expected improve- 
ments in the accuracy in the Earth-Mercury ranging. At present, the 1-way 
range residuals of Mercury from radar-ranging span 30 yr (1967-1997) and are 
at a few km- level (Figure B-2 of Folkner et al. (2008)); the same holds for the 
1-way Mercury radar closure residuals covering 8 yr (1989-1997, Figure B-3 a) 
of Folkner et al. (2008)). There are also a pair of Mariner 10 range residuals 
in the 70s at Mercury at 0.2 km level (Figure B-3 b) of Folkner et al. (2008)). 

Here we outline the general approach followed. In order to numerically ob- 
tain the effect of a given gravitational perturbation P of the Newtonian Sun's 
monopole on the range |p| between the Earth-Moon Barycenter (EMB) and 
a given planet, we used MATHEMATICA to simultaneously integrate with 
the Runge-Kutta method the equations of motion in Cartesian coordinates 
of EMB and the planet considered with and without the perturbation P by 
using the same set of initial conditions. We adopted the ICRF/J2000.0 ref- 
erence frame, with the ecliptic and mean equinox of the reference epoch, i.e. 
J2000, centered at the Solar System Barycenter (SSB); the initial conditions at 
the epoch J2000 were retrieved with the HORIZONS WEB interface by JPL, 
NASA. The temporal interval of the numerical integration has been taken 
equal to At = 2 yr because most of the present-day available time series of 
range residuals approximately cover similar temporal extensions; moreover, 
also the typical operational time spans forecasted for future PLR technique 
are similar. The basic model adopted consists of the barycentric equations of 
motion of the Sun, the eight planets, the Moon, Ceres, Pallas, Vesta, Pluto 
and Eris, to be simultaneously integrated; the forces acting on them include 
the mutual Newtonian N— body interactions, the perturbation due to the so- 
lar quadrupolar mass moment J2, the effect of two rings modeling the actions 
of the minor asteroids and of the Trans- Neptunian Objects (TNOs), and the 
post-Newtonian GEM field of the Sun (with its general relativistic value). 

In order to preliminarily evaluate the potential measurability of the effects 



3 Its cheaper version ASTROD I makes use of one spacecraft in a Venus-gravity- 
assisted solar orbit, ranging optically with ground stations (Appourchaux et al., 
2009). 
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considered, the computed differences A|p| = |pp| — \pr\, where R refers to a 
reference orbit which does not contain the perturbation P, were subsequently 
compared to the available time series of the range residuals for the inner 
planets which set the present-day accuracy level in ranging to planets (Folkner 
et al., 2008; Pitjeva, 2010). About the possibility that a given, unmodeled 
dynamical effect may show or not its signature in the range residuals, it must 
be considered that the magnitude of such an effect should roughly be one order 
of magnitude larger than the range residuals accuracy. This to avoid the risk 
that it may be absorbed and partially or totally removed from the signature 
in the process of estimation of the initial conditions and of the other numerous 
solve-for parameters in the real data reduction. 

Depending on the dynamical effect one is interested in, some of the per- 
turbations examined here are to be considered as sources of noise inducing 
systematic bias on the target signal. For example, if the goal of the analysis 
is, as in our case, the LT effect, then the range perturbation due to the TNOs 
is clearly a source of potential systematic error which has to be evaluated. 
Thus, our plots are useful to assess the level of aliasing of several potential 
sources of aliasing for some non-Newtonian effects and the correlations that 
may occur in estimating them. Dynamical effects which are viewed as noise in 
a given context can also be regarded as main targets in another one; see, e.g., 
the proposed determination of asteroid masses through the ASTROD mission 
(Su et al., 1999). In this respect, the LT signature, if not properly modeled, 
should be regarded as a source of potential bias. 



2. 1 The gravitomagnetic, Lense- Thirring field of the Sun 

Figure 1 depicts the Earth-Mercury range perturbation due to the Sun's 
GM field, neither considered so far in the dynamical force models of the plan- 
etary ephemerides nor in the BepiColombo analyses. It is important to note 
that the value of eq. (5) for the Sun's angular momentum does not come from 
planetary orbital dynamics, so that there is no risk of a-priori "imprinting" of 
GTR itself on range tests of the solar LT effect which could, thus, be regarded 
as genuine and unbiased. 

As a small technical note, we mention that we rotated the reference frame 
to the mean ecliptic at the epoch to the Sun's equator by the Carrington angle 
% = 7.15 deg (Beck & Giles, 2005) because eq. (3) holds in a frame with its z 
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axis aligned with S. 

The peak-to-peak amplitude of the LT signal is up to 17.5 m over 2 yr, 
which, if on the one hand is unmeasurable from currently available radar- 
ranging to Mercury, on the other hand corresponds to a potential relative 



7 



EMB -Mercury range: SSB numerical calculation 
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Fig. 1. Difference A\p\ = \pp\ — \p^\ in the numericalfy integrated EMB-Mercury 
ranges with and without the perturbation due to the Sun's GM field over At = 2 
yr. The same initial conditions (J2000) have been used for both the integrations. 
The state vectors at the reference epoch have been retrieved from the NASA JPL 
Horizons system. The integrations have been performed in the ICRF/J2000.0 ref- 
erence frame, with the reference {xy} plane rotated from the mean ecliptic of the 
epoch to the Sun's equator, centered at the Solar System Barycenter (SSB). 

accuracy in measuring it with BepiColombo of 2 — 5 x 10~ 3 ; this clearly shows 
that the solar GM field should be taken into account in future analyses and 
data processing. Otherwise, it may alias the recovery of other effects, as it will 
become clearer later. 



2.2 The gravitoelectric, Schwarzs child field of the Sun 

In Figure 2 we plot the effect of the GE field of the Sun on the Earth- 
Mercury range. 

Figure 2 can be compared with Figure 1 of Milani et al. (2010), obtained 
for unspecified initial conditions 4 : they are quite similar. 

The maximum variation of the signal is of the order of 4 x 10 5 m, cor- 
responding to a measurement accuracy of about 2.5 x 10~ 7 . The expected 
realistic accuracy in determining j3 and 7 is 2 x 10 -6 in BepiColombo (Milani 
et al., 2002). To this aim, let us note that, since the LT effect depends on 7, 
neglecting it may alias the determination of 7 through the larger GE signal 
at 4 x 10" 5 level. 



4 It also includes the Shapiro delay contribution. 
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EMB -Mercury range: SSB numerical calculation 
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Fig. 2. Difference A|p| = \pp\ — \p&\ in the numericalfy integrated EMB-Mercury 
ranges with and without the perturbation due to the Sun's Schwarzschifd fiefd over 
At = 2 yr. The same initial conditions (J2000) have been used for both the integra- 
tions. The state vectors at the reference epoch have been retrieved from the NASA 
JPL Horizons system. The integrations have been performed in the ICRF/J2000.0 
reference frame, with the ecliptic and mean equinox of the reference epoch, centered 
at the Solar System Barycenter (SSB). 

2.3 The Newtonian effect of the oblateness of the Sun 



Figure 3 shows the nominal effect of the Sun's quadrupolar mass moment 
on the Earth-Mercury range for J 2 = 2 x 10~ 7 . Its action has been modeled 
as (Vrbik, 2005) 



.4 



■)■> 



3J 2 R 2 GM 
2V A 



1 — 5 (f • k 



f + 2 ( f • k ) k 



where R is the Sun's mean equatorial radius and k is the unit vector of the z 
axis directed along the body's rotation axis. 

Since eq. (8) holds in a frame with its {xy} plane coinciding with the body's 
equator, also in this case we rotated the the mean ecliptic at the epoch to the 
Sun's equator. 

The signal of Figure 3 has a maximum span of 300 m, corresponding to an 
accuracy measurement of 3 x 10~ 4 . A determination of the solar J 2 accurate 
to 10~ 2 is one of the goals of BepiColombo (Milani et al., 2002); knowing 
precisely J 2 would yield important insights on the internal rotation of the 
Sun. At present, the uncertainty in it is about 10% (Fienga et al., 2010). The 
solar quadrupole mass moment may play the role of source of systematic bias 
with respect to, e.g., some non- Newtonian dynamical effects. Concerning the 
GE signal previously analyzed, the mismodeled J 2 signature would impact it 
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EMB-Mercury range: SSB numerical calculation 
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Fig. 3. Difference A|p| = \pp\ — in the numerically integrated EMB-Mercury 
ranges with and without the nominal perturbation due to the Sun's quadrupole 
mass moment J2 = 2.0 x 1CP 7 over At = 2 yr. The same initial conditions (J2000) 
have been used for both the integrations. The state vectors at the reference epoch 
have been retrieved from the NASA JPL Horizons system. The integrations have 
been performed in the ICRF/J2000.0 reference frame, with the reference {xy} plane 
rotated from the mean ecliptic of the epoch to the Sun's equator, centered at the 
Solar System Barycenter (SSB). 

a 7.5 x 10 -5 level. It is important to note that the patterns of the two signals 
are rather different. Conversely, the determination of J 2 at the desired level of 
accuracy may be affected by other unmodeled/mismdeled dynamical effects 
acting as systematic sources of aliasing on it. For example, the GM signature 
may affect the determination of J2 at 12% level. On the other hand, in order 
to allow for a determination of the LT effect itself, the Sun's quadrupole mass 
moment should be known with an accuracy better than the present-day one by 
at least one order of magnitude; this is just one of the goals of BepiColombo. 
Anyway, also the GM and the J 2 patterns are different. 



2.4 The Newtonian effect of the ring of the minor asteroids and Ceres, Pallas 
and Vesta 



In Figure 4 we depict one potential source of systematic bias, i.e. the action 
of the ring of minor asteroids (Fienga et al., 2010). We modeled it following 
Fienga et al. (2008). In general, for those planets for which r < -R r i ng , by 
posing a = r/i? ring , we obtained 

outer rmg " 2rm mp { + 8 64 ) ' K } 
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EMB -Mercury range: SSB numerical calculation 
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Fig. 4. Difference A|p| = \pp\ — \pn\ in the numerically integrated EMB-Mercury 
ranges with and without the nominal perturbation due to the ring of minor asteroids 
with m r i ng = 1 x lCT 10 Mo (Fienga et al., 2010) and i? r i ng = 3.14 au (Fienga et al., 
2010) over At = 2 yr. The same initial conditions (J2000) have been used for 
both the integrations. The state vectors at the reference epoch have been retrieved 
from the NASA JPL Horizons system. The integrations have been performed in the 
ICRF/J2000.0 reference frame, with the ecliptic and mean equinox of the reference 
epoch, centered at the Solar System Barycenter (SSB). 
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(10) 



Recall that the Laplace coefficients are defined as 

&W( a ) = I / cos ^ (11) 

7r J (1 — 2acosifj + a 2 ) s 



o 



where s is a half-integer; a useful approximate expression in terms of a series 
can be found in Murray & Dermott (1999, p. 237) 



U(j) _ s(s+l)...(s+j-l) j 
u s — l-3...j " 



1 i s(s+j) 2. 
1+ (1+3) a 



, s(s+l)(s+j)(s+j+l) 4 
^ l-2(i+l)(i+2) « 



By assuming for the ring of the minor asteroids a nominal mass of m ring = 
1 x 1O~ 1O M (Fienga et al., 2010) and a radius -R ring = 3.14 au (Fienga 
et al., 2010), it would impact the Mercury range at 4 m level (peak-to- 
peak amplitude), which will be, in fact, measurable. Its nominal bias on the 
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EMB -Mercury range: SSB numerical calculation 
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Fig. 5. Difference A|p| = \pp\ — \pp\ in the numericaliy integrated EMB-Mercury 
ranges with and without the nominai perturbation due to Ceres, Pallas, Vesta (Pit- 
jeva & Standish, 2009) over At = 2 yr. The same initial conditions (J2000) have 
been used for both the integrations. The state vectors at the reference epoch have 
been retrieved from the NASA JPL Horizons system. The integrations have been 
performed in the ICRF/J2000.0 reference frame, with the ecliptic and mean equinox 
of the reference epoch, centered at the Solar System Barycenter (SSB). 



Schwarzschild, J 2 and LT signals would be 1 x 10~ 5 , 1.3 x 10~ 2 ,2.3 x KT 1 , 
respectively. Anyway, the present-day level of uncertainty in the mass of the 
ring is <5m ring = 0.3 x 1O~ 1O M (Fienga et al., 2010). Thus, the impact of such 
a mismodeling would be, 3 x 10~ 6 ,4 x 10~ 3 , 7 x 10~ 2 , respectively; it cannot 
be considered negligible. 



The effect of Ceres, Pallas and Vesta on the determination of some Newto- 
nian and non- Newtonian parameters with BepiColombo has been preliminarily 
investigated in Ashby et al. (2007). Here in Figure 5 we show the nominal per- 
turbation on the Earth-Mercury range due to the combined actions of Ceres, 
Pallas and Vesta; the values for their masses have been retrieved from Pitjeva 
& Standish (2009). Its peak-to-peak amplitude amounts to 80 m; thus, their 
signature would be measurable at a 0.6 — 1 x 10 -3 level. Anyway, the mismod- 
eled solar quadrupole mass moment would bias their signal at 4 x 10 _1 level. 
The LT effect, if unmodeled, would have an impact at 2.2 x lO^ 1 level. The 
present-day relative uncertainties in their masses are 6 x 10~ 3 , 3 x 10~ 2 , 2 x 10~ 2 
respectively (Pitjeva & Standish, 2009). This implies a mismodeled signal with 
a peak-to-peak amplitude of 50 cm. It would impact the Schwarzschild, J 2 and 
LT range perturbations at 1 x 10~ 6 , 2 x 10 -3 , 3 x 10~ 2 level, respectively. 
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EMB -Mercury range: SSB numerical calculation 
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Fig. 6. Difference A\p\ = \pp\ — \p^\ in the numericalfy integrated EMB-Mercury 
ranges with and without the nominaf perturbation due to the ring of Trans-Nep- 
tunian Objects with m r i ng = 5.26 x 10~ 8 M Q (Pitjeva, 2010) and i? r i ng = 43 au 
(Pitjeva, 2010) over At = 2 yr. The same initial conditions (J2000) have been 
used for both the integrations. The state vectors at the reference epoch have been 
retrieved from the NASA JPL Horizons system. The integrations have been per- 
formed in the ICRF/J2000.0 reference frame, with the ecliptic and mean equinox 
of the reference epoch, centered at the Solar System Barycenter (SSB). 

2. 5 The Newtonian effect of the Trans- Neptunian Objects 

The situation is different for another potential source of systematic un- 
certainty, i.e. the Trans- Neptunian Objects (TNOs). Figure 6, obtained by 
modeling them as a ring with m ring = 5.26 x 1O~ 8 M (Pitjeva, 2010) and 
-Rrmg = 43 au (Pitjeva, 2010), shows that their maximum effect would amount 
to 80 cm. We used the same formulas as for the asteroid ring. Such an effect, 
not taken into account so far, would be better measurable than that by the 
minor asteroids. This implies a bias of 2 x 10 -6 on the Schwarzschild signal, 
3 x 10 -3 for J 2 and 4.5 x 10~ 2 for the Lense-Thirring effect. A major concern 
is that the mass of the TNOs is far from being accurately known, so that an 
uncertainty as large as 100% should be applied. 



3 Summary and Conclusions 

In view of a possible future implementation of some interplanetary laser 
ranging facilities accurate to cm-level (4.5 — 10 cm for BepiColombo in 1 — 8 
yr), we have numerically investigated how the ranges between the Earth and 
Mercury are affected by certain Newtonian and non-Newtonian dynamical 
effects by simultaneously integrating the equations of motion of all the major 
bodies of the Solar System plus some minor bodies of it (Ceres, Pallas, Vesta) 
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in the SSB reference frame over a time span two years long. 

It turns out that the general relativistic gravitomagnetic Lense-Thirring 
effect of the Sun, not modeled so far either in the planetary ephemerides or in 
the analyses of some spacecraft-based future missions like, e.g., BepiColombo, 
does actually fall within the measurability domain of future cm-level ranging 
facilities. The more favorable situation occurs for Mercury because the relative 
measurement accuracy is of the order of 2 — 5 x 10~ 3 by assuming a 4.5 — 10 
cm uncertainty in the Earth-Mercury ranging, as expected for BepiColombo 
over some years of operations. 

If not properly modeled and solved-for, the Lense-Thirring effect may also 
impact the determination of other Newtonian and post-Newtonian parame- 
ters at a non-negligible level, given the high accuracy with which their mea- 
surement is pursued. For example, in the case of BepiColombo the expected 
accuracy in determining 7 and (5 from the range perturbation due to the 
Schwarzschild field of the Sun is of the order of 10~ 6 ; the Lense-Thirring range 
perturbation would impact the Schwarzschild one at 4 x 10 -5 level. Another 
goal of the BepiColombo mission is a measurement of the Sun's quadrupole 
mass moment J 2 accurate to 10~ 2 ; the unmodeled Lense-Thirring effect would 
bias it at 10 _1 level. 

From the point of view of a measurement of the Sun's gravitomagnetic field 
itself, it results that a major concern would be the solar oblateness; it should 
be known at a 1CT 2 level of accuracy-which is just the goal of BepiColombo-to 
allow for a reduction of its aliasing impact on the Lense-Thirring signal down 
to just 17%. The ring of the minor asteroids should be taken into account 
as well because its mismodeling would impact the gravitomagnetic signal at 
about 7 x 10 -2 . The lin gering uncertainty in the masses of Ceres, Pallas, Vesta 
translates into a potential bias of about 4.5 x 10~ 2 — 1CT 3 . The TNOs, not 
modeled so far apart from the EPM ephemerides, would nominally affect it 
at a 4.5 x 10~ 2 level; it must be considered that there is currently a high 
uncertainty in their mass. However, it must be noted that the previous figures 
have been obtained by comparing the peak-to peak amplitudes of the various 
time-dependent range signals. Actually, the time signatures of such sources 
of systematic bias are different with respect to the gravitomagnetic one; this 
would greatly help in disentangling it from the noisy effects. 

Table 2 summarizes our findings. 

Finally, we point out that the analysis presented here should be regarded 
just as a preliminary investigation. Actually, extensive simulations should be 
performed in which a gravitomagnetic-dedicated parameter should be esti- 
mated along with the wealth of other ones which are routinely solved-for in 
the data reduction procedures in order to check how the relativistic signal of 
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interest may be affected by the estimation of, say, the initial conditions of 
the planets/spacecraft whose data are processed. Indeed, it must be recalled 
that an unmodelled dynamical effect may be, partially or totally, absorbed 
in some of the estimated parameters, especially if its magnitude is not large 
enough with respect to the characteristic accuracy level of the observations. 
Anyway, the implementation of such a non-trivial task is beyond the scopes of 
the present paper, and may well constitute the subject for further researches. 
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Table 1 

First line: uncertainties (in m) in the average heliocentric distances of the inner 
planets obtained by propagating the formal errors in a and e according to Table 3 
of Pitjeva (2008); the EPM2006 ephemerides were used by Pitjeva (2008). Second 
line: maximum differences (in m) between the EPM2006 and the DE414 (Standish, 
2006) ephemerides for the inner planets in the time interval 1960-2020 according 
to Table 5 of Pitjeva (2010). Third line: maximum differences (in m) between the 
EPM2008 (Pitjeva, 2010) and the DE421 (Folkner et al., 2008) ephemerides for the 
inner planets in the time interval 1950-2050 according to Table 5 of Pitjeva (2010). 
They have to be compared with the characteristic gravitomagnetic length of the 



Sun If = 319 m. 



Type of orbit uncertainty 


Mercury 


Venus 


Earth 


Mars 


5 (r) (EPM2006) 


38 


3 


1 


2 


EPM2006-DE414 


256 


131 


17.2 


78.7 


EPM2008-DE421 


185 


4.6 


11.9 


233 



Table 2 

Maximum peak-to-peak nominal amplitudes, in m, of the Earth-Mercury range sig- 
nals over At = 2 yr due to the dynamical effects listed. We adopted the standard 
value = 2.0 x 10~ 7 (Fienga et al., 2010) for the quadrupole mass moment of the 
Sun. It is presently known at a 10% level of accuracy. For its proper angular momen- 
tum we used Sq = 190.0 x 10 39 kg m 2 s _1 (Pijpers, 1998, 2003) from helioseismology. 
For the ring of the minor asteroids we used m r i Ilg = (1±0.3) x 10~ 10 Mq, -R r i ng = 3.14 
au (Fienga et al., 2010), while for the TNOs, modeled as massive ring as well, we 
adopted m ring = 5.26 x 1O- 8 M , i^ng = 43 au (Pitjeva, 2010). The masses of 
the major asteroids Ceres Pallas, Vesta, accurate to 10~ 2 — 10 -3 level, have been 
retrieved from Pitjeva & Standish (2009). 

Dynamical effect Peak-to-peak amplitude (m) 

Solar Schwarzschild 4 x 10 5 

Solar J® 300 

Ceres, Pallas, Vesta 80 

Solar Lense-Thirring 17.5 

Ring of minor Asteroids 4 

TNOs 0.8 
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